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Abstract 

The properties of triaxial bright solitons (TBSs) made of attractive Bose- 
Einstein condensed atoms under transverse anisotropic harmonic confine- 
ment are investigated by using a variational approach. We show that these 
metastable TBSs change their shape from cigar-like to disc-like by increasing 
the inter-atomic strength. Moreover, we find that the collective oscillations 
of a TBS strongly depend on the anisotropy parameter of the external po- 
tential. We calculate in detail the properties of TBS close to the collapse, 
which sets up at a critical value of the inter-atomic strength. This critical 
strength is maximal in the isotropic (axially symmetric) case and slightly re- 
duces by increasing the anisotropy. Finally, we investigate the formation of 
multiple TBSs via modulational instability induced by a sudden change of the 
scattering length from positive to negative. 
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Bright solitons (BSs) have been recently obtained in two experiments with Bose-Einstein 
condensates (BECs) of 7 Li dilute vapors at ultra-low temperatures [1,2]. In these experi- 
ments the attractive BEC is confined in two directions by a transverse cylindric isotropic 
harmonic potential and travels in the third direction, the longitudinal cylindric axis, without 
a relevant spreading. The theoretical analysis of these BS configurations has been devel- 
oped by various authors [3-8]. It has been shown that the confinement in two directions is 
necessary to create single [3-6] or multiple [7,8] metastable BSs, which collapse above a crit- 
ical number of particles [5,6,8]. In the previous theoretical investigations only an isotropic 
transverse cylindric harmonic confinement has been studied. Under these conditions the BS 
is cigar-shaped, it becomes less cigar-shaped by increasing the inter-atomic strength and it 
acquires a quasi-spherical shape only very close to the critical inter-atomic strength of the 
collapse [5,6]. 

In this paper we consider the case of an attractive BEC under anisotropic harmonic 
confinement. This condensate BS is triaxial and we find that it can be cigar-shaped or disc- 
shaped by changing the inter-atomic strength. The triaxial bright soliton (TBS) collapses 
at a critical inter-atomic strength that decreases by increasing the anisotropy parameter 
of the transverse harmonic potential. By using a variational approach we calculate the 
collective oscillations of the TBS. We find that close to the collapse the longitudinal collective 
frequency reaches its maximum value and then quickly goes to zero. The maximum value 
of the longitudinal collective frequency strongly depends on the anisotropy parameter. We 
calculate also the number of TBSs which can be obtained, starting with a repulsive and 
axially uniform BEC, by suddenly changing the scattering length from positive to negative. 
This mechanism, known as modulational instability [7,8], has been experimentally applied 
to get a soliton train of axially symmetric BSs [2]. 

Nowadays BECs are routinely created and trapped in magnetic or optical traps. Such 
traps can be modeled with good accuracy with harmonic potentials. Here we study an 
attractive BEC with N atoms confined in the cylindric transverse direction by an anisotropic 
harmonic potential U(t) = {m/2){u)\x 2 + u>%y 2 ), where m is the mass of a Bose-condensed 
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atom and u\ and c^2 are the two frequencies of the transverse confinement. The ratio 
A = uj 2 /uji is the anisotropy parameter of the harmonic potential. Note that there is no 
confinement along the z axis. 

The scaled Lagrangian density of this attractive BEC is given by 



where ip(r,t) is the macroscopic wave function of the BEC normalized to N, a s is the s- 
wave scattering length (a s < with attractive atoms) and an = ih/imajn)) 1 ^ 2 is the average 
harmonic length of the trapping potential with uo H = (c^ic^) 1 / 2 - Note that in Eq. (1) lengths 
are in units of a#, energies are in units of Hcuh, and the time is in units of ujJj 1 . In the rest 
of the paper we use these scaled variables. 

The Euler-Lagrange equation obtained from the Lagrangian density C is the familiar 
3D time- dependent Gross-Pitaevskii equation (GPE), which describes very accurately a 
pure BEC of Bosonic vapors [9]. By using the 3D GPE one finds that the ground-state of 
an attractive BEC is a collapsed state of zero radius and infinite negative energy but, as 
previously stressed, metastable states can be obtained, up to a critical number of particles, 
under external confinement in two or three directions [5,6]. The stationary metastable states 
of an attractive BEC under confinement in two directions are referred as BSs because they 
are self-confined in the third direction and can propagate without spreading along that 
direction [10]. For A = 1 one has axial BSs, whose statical and dynamical properties have 
been analyzed in detail in previous papers [5-8]. For A ^ 1 one gets TBSs, whose properties 
are investigated in this paper. 

Instead of numerically solving the full 3D GPE we investigate the properties of the TBS 
by using a variational approach [11], which is computationally much less demanding. Our 
variational ansatz for the macroscopic wave function of the attractive BEC is the following 



with r = {x\,X2,xz) = (x,y,z). <Jk(t) and ftkit) are the time-dependent variational parame- 
ters. The a kit) are the widths of the TBS in the three axial directions. Note that in order 




(1) 
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to describe the time evolution of the variational wave function, the phase factor i(5k{t)x\ is 
needed. 

By using this Gaussian trial wave function, after spatial integration of the Lagrangian 
density C of Eq. (1), one finds a new effective Lagrangian L that has <J k (t) and f3k(t) as 
generalized coordinates. The three Euler-Lagrange equations of L with respect to a k (t) are 
given by 

^2 l/2 7 

af V 7T a{ a 2 a 3 
o- 2 + A 2 a 2 = — - y 2 — ' ( 3 ) 

(7| V 7T <7i<7|<7 3 
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where 7 = 7V|a s |/aff is the scaled inter-atomic strength and Ai = oo\/ooh = 1/a/A and 
A2 = = \/A are the scaled frequencies of the transverse harmonic potential. The three 

Euler-Lagrange equations of L with respect to /9fc(t) give instead the following expressions: 
Pk — —&kl (2<7 fc ), with k = 1,2, 3. Thus the time dependence of /^(t) is fully determined by 
that of (Jfc(t). 

The stationary metastable states of Eqs. (3), namely the TBS static configurations, are 
found by setting a k = 0. The resulting algebraic equations can be easily solved numeri- 
cally. The solutions correspond to metastable states if they are local minima of the effective 
potential energy 

+ i (X\a\ + \y 2 ) - ^--L- (4) 

associated to the Newtonian Eqs. (3). In Figure 1 we plot the three widths <j\, a 2 and 
<T 3 of the metastable TBS as a function of the inter-atomic strength 7. Note that the 
transverse widths <J\ and a 2 do not change very much with respect to the non-interacting 
values o\ = \/y/\[ = A 1 / 4 and a 2 = l/\ /r X 2 = 1/A 1 / 4 , while the longitudinal width <t 3 is 
divergent for 7 = and approaches one of the two transverse widths for 7 close to the critical 



value 7c. For 7 > 7 C there are no more metastable solutions: one has the so-called collapse 
of the attractive BEC. 

It is important to observe that, by using the Gaussian variational approach, the critical 
value 7 C of the collapse is slightly overestimated with respect to the "exact" numerical 
calculations. For instance, with A = 1 the "exact" result gives 7 C = 0.676 [12] while our 
variational method predicts 7 C = 0.778. Moreover, in the case of asymmetric harmonic trap 
along the three axes, it has been shown that the relative error of 7 C by using the variational 
approach is always within 10% [12]. 

Figure 1 shows that for small values of 7 the TBS is highly cigar-shaped but near 7 C the 
TBS becomes disk-shaped due to the enormous compression along the longitudinal z axis. 
In particular, close to the collapse with A = 1 the system is spherical-shaped but choosing 
a large anisotropy (A ^> 1 or A 1) it is strongly disk-shaped. Note that, due to the 
symmetry of the problem, one has o"i(A) = o" 2 (A _1 ) and (73(A) = o" 3 (A _1 ). 

The Gaussian approximation of the TBS wave function can be used to study also the 
collective oscillations the TBS. The diagonalization of the Hessian matrix d 2 W / da^dai of 
the effective potential energy W of Eq. (4) gives three frequencies f^, f2 2 and f2 3 of col- 
lective excitations around the TBS solution. The 3x3 Hessian matrix can be numerically 
diagonalized by choosing, for fixed 7 and A, the widths that satisfy the Eqs. (3) with 
a k = 0. 

In Figure 2 we plot such frequencies. Only for 7 = the frequencies fli, fl 2 and fl 3 can be 
interpreted as collective oscillations along x axis, y axis and z axis, respectively; nevertheless 
the mixing angle remains quite small also for finite values of 7 so they can be associated 
to the motion along the three axes. The transverse frequencies fii and f2 2 are practically 
constant with respect to 7 and equal to the non-interacting values Qi = 2Ai = 2/y/X and 
f2 2 = 2A 2 = 2\f\; only very close to 7 C the two transverse collective frequencies suddenly 
grows. Instead, the longitudinal frequency f2 3 is zero for 7 = 0, it increases with 7 but close 
to 7c it falls down to zero. Note that, again due to the symmetry of the problem, one has 
Sli(A) = fi 2 (A~ 1 ) and S1 3 (A) = S7 3 (A _1 ). The maximum value fig" 2 -* of the frequency f2 3 is 
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plotted in Figure 3 (dot-dashed line) as a function of A. In Figure 3 is also shown (full line) 
the graph of 7 C (A). Both and 7 C slowly decreases by increasing the anisotropy (note 
the logarithmic scale in the A axis of Figure 3). 

In a recent experiment [2] a train of bright solitons has been obtained with an initially 
attractive and axially uniform BEC by a sudden change of the scattering length a s from 
positive (repulsive) to negative (attractive), induced by an external magnetic field (Feshbach 
resonance [13]). The formation of this multi-soliton configuration can be explained as due 
to the modulational instability of the time-dependent wave function of the BEC, driven by 
imaginary Bogoliubov excitations [7,8]. An analytical formula for the number N s of bright 
solitons generated with a quasi-lD condensate via modulational instability has been derived 
by Al Khawaja et al. [7] and by Salasnich, Parola and Reatto [8]. The formula has been 
generalized by Salasnich [8] to the case of a 3D BEC under isotropic transverse confinement. 
Here we predict that the anisotropy of the transverse harmonic confinement strongly affects 
the number N s of bright solitons. 

Let us consider a BEC made of N atoms with positive scattering length (a s > 0) trapped 
by a harmonic potential in the transverse direction (x, y) and by a box potential of length L 
in the longitudinal axial direction (z). In this way the stationary BEC is uniform along the 
z axis (axially uniform) and it can be modeled by the following variational wave function 



where n = N/ L is the axial density of the BEC and the two variational parameters U\ and 
<t 2 are the width of the BEC along x and y axes. It is straightforward to derive the energy 
of this BEC from the Lagrangian density of Eq. (1) by setting d t ip = 0, inserting the Eq. 
(5) into Eq. (1) and integrating the lagrangian density over x and y. In this way one finds 
for the energy the following formula 



where g = 2a s /an- This energy depends on the two variational parameters o\ and o<i- The 
minimization of the energy E with respect to <j\ and <j 2 gives the equations 
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Ai<7i = 1 + #n— , A 2 <t 2 = 1 + gn— . 



(7) 




From these nontrivial equations one numerically finds, for a fixed inter- atomic strength g, the 
transverse widths a k of the Bose gas as a function of the axial density n. Having derived the 
functions Ofc(n), one can numerically determine the chemical potential \i = dE/dn. From 
the chemical potential ji one gets the sound velocity c in the longitudinal axial direction, 



the energy E, the chemical potential ji and the sound velocity c obtained with Eqs. (6-7) 
give precisely the ID GPE results in the weak-coupling ID limit (a± = A 1 / 4 and a 2 = 1/A 1 / 4 ) 
while they have the correct (Thomas-Fermi) 3D GPE behaviour in the strong-coupling 3D 
limit (a 4 = \ 2 gn and a\ = gn/\ 2 ) (see also [14] for the A = 1 case). 

The sound velocity is useful to calculate the Bogoliubov excitations e k = 
[(k 2 /2) (k 2 /2 + 2c 2 )] 1 / 2 of the axially uniform BEC. By suddenly changing the scattering 
length a s to a negative value, the excitations frequencies corresponding to k < k c = 2\c\ be- 
come imaginary and, as a result, small perturbations grow exponentially in time. The maxi- 
mum rate of growth is at k = k c j y/2 and the wavelength of this mode is A = 2n/k . The ra- 
tio L/\ gives an estimate of the number N s of TBSs which are generated: N s = \c\L/ (y/2ir), 
where \c\ is now a function of n\g\ = with 7 = N\a s \/au- Note that the criterion 

of modulational instability we use has been recently re-derived with a sophisticated time- 
dependent analysis [15]. 

In Figure 4 we plot N s /L as a function of 7/L for different values of the anisotropy 
parameter A. For small values of 7/L we are in the weak-coupling quasi-lD regime but 
approaching 7/L = 1/2 our procedure gives a divergent number of solitons because one 
of the two transverse widths shrinks to zero for a large attractive inter-atomic interaction. 
Figure 4 shows that the number N s of TBSs generated via modulational instability increases 
with the anisotropy parameter A. Moreover, as expected, we find N S (X) = iV^A -1 ). 

In conclusion, we have analyzed metastable states (triaxial bright solitons) and collec- 
tive oscillations of an attractive Bose condensate under anisotropic harmonic confinement. 



which satisfies the equation c = (nd/i/dn) 



x l 2 = [nd 2 E / dn 2 ) 1 / 2 . It is important to stress that 
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We have predicted for the triaxial bright soliton a transition from a cigar-like shape to a 
disk-like shape by increasing the inter-atomic strength, namely the number of particles or 
the (negative) scattering length, for instance by using Feshbach resonances. We have also 
shown that the number N s of triaxial bright solitons generated with a sudden change of 
the scattering length from positive to negative grows with the anisotropy of the transverse 
confinement. Moreover N s diverges at a critical inter- atomic strength, where one of the two 
transverse widths of the soliton train shrinks to zero. 
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FIG. 1: Widths of the TBS as a function of the inter-atomic strength 7 = N\a s \/an, where 
A = U2/0J1 = A2/A1 and an = (^/{itlujh)) 1 ^ 2 with ujh = (wi^) 1 ^ 2 - Dotted line: <j\, dashed 
line: 02, solid line: 03. Lengths are in units of the harmonic length an of the external transverse 
potential. 
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FIG. 2: Collective frequencies of the TBS as a function of the inter-atomic strength 7 = 
N\a s \/a,H for some values of the anisotropy parameter A = U2/0J1 = A2/A1. Frequencies are in 
units of the harmonic frequency ujh = (^l^) 1 ^ 2 of the external transverse potential. 
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FIG. 3: Properties of the TBS at the collapse as a function of the anisotropy parameter 
A = W2/0J1 = A2/A1. The solid line is the critical inter-atomic strength j c . The dot-dashed line is 
the maximum value of the collective frequency Q3 along the longitudinal z axis. Units as in 
Fig. 2. 



13 




0.1 0.2 0.3 0.4 0.5 

y/L 

FIG. 4: Number N s of bright solitons generated via modulational instability as a function of 
the inter-atomic strength 7 = N\a s \/a,H- L is the initial axial length of the BEC with N atoms. 
\ = UJ2/0J1 = A2/A1 is the anisotropy parameter of the transverse harmonic confinement. Units as 
in Fig. 1. 
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